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ABSTRACT 


Data  from  the  short-period  vertical  subarray  near  Oyer  and  the  multicomponent 
long-period  seismograph  near  Faldalen,  Oyer  and  Trysil,  Norway,  were  used  to  in¬ 
vestigate  the  discrimination  of  earthquakes  and  underground  nuclear  detonations  using 
long-period  to  short -period  energy  ratio  (M^  v£  m^). 


iii 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 


INTRODUCTION 


Data  from  the  short-period  vertical  subarray  near  Oyer  and  the  multicomponent 
long-period  seismographs  near  Faldalen,  Oyer  and  Trysil,  Norway,  were  used  to  in¬ 
vestigate  the  discrimination  of  earthquakes  and  underground  nuclear  detonations  using 
long -period  to  short -period  energy  ratio  (M  vs  m, ).  These  stations  were  installed 
as  part  of  the  site  survey  for  the  Norwegian  Seismic  Array  (NORSAR). 

To  insure  first  that  the  long-period  instruments  were  operating  normally  and 
that  the  data  is  comparable  with  that  recorded  elsewhere,  the  surface  wave  magnitude 

(M  )  recorded  at  LASA  and  Faldalen  was  studied  for  a  population  of  21  events,  and 
s 

then  a  further  comparison  was  made  between  the  three  Norway  sites  using  a  population 
recorded  at  at  least  two  of  them. 

The  Mg  vs  m^  relation  for  four  populations  of  events  recorded  at  LASA  was 
compared  with  initial  data  from  Norway.  The  surface  wave  magnitude*  was  computed 
from  the  amplitude  of  the  20-second  period  Rayleigh  wave  using 

M  =  log  A  —  log  B 
s 

where  A  is  the  zero-to-peak  amplitude  of  the  20-second  Rayleigh  wave,  in  millimicrons, 
recorded  on  a  long-period  vertical  instrument  and  log  B  is  determined  from  epicentral 
distance  by 


log  B  =  1.182  -1.656  log  (A) 
where  A  is  the  epicentral  distance  in  degrees. 
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The  body-wave  magnitude  (m^)  was  determined  from  the  amplitude  of  the 

2 

short -period  vertical  wave  recorded  at  teleseismic  distances,  using 

mb  =  log  W/T  +  Q  +  S 

where  W  is  the  maximum  zero -to -peak  amplitude  (in  millimicrons)  of  the  first  four 

cycles  of  the  P-wave,  T  is  the  period  (in  seconds)  of  the  short-period  cycle  used,  Q 

depends  on  epicentral  distance  and  is  tabulated  in  Fig.  5  of  a  paper  by  Gutenberg  and 
2 

Richter,  and  S  is  a  station  correction  which  is  taken  to  be  zero  for  convenience. 

With  the  commencement  of  data  acquisition  in  Norway  we  had  an  opportunity  to 
investigate  the  western  United  States  region  at  teleseismic  distances.  This  permitted 
a  study  of  the  character  of  Mg  vs.  m^,  as  well  as  its  effectiveness  as  a  discriminant, 
for  events  from  this  region.  The  data  obtained  provides  grounds  for  some  speculation 
on  the  peculiarities  of  the  western  United  States. 

3 

The  Mg  and  m^  values  were  computed  on  the  data  analysis  console.  This 
approach  offered  precise  determination  of  Mg  and  while  permitting  maximum 
human  interaction  when  analyzing  the  data,  e.  g.  bad  channels  could  be  detected,  and 
time  delays  for  beams  could  be  theoretically  determined  or  hand  picked. 
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INTER-SITE  COMPARISON  OF  SURFACE  WAVE  MAGNITUDE 


Twenty-one  events  were  found  with  surface  waves  observed  at  both  Faldalen 

and  LASA.  Figure  1  shows  the  M  at  NORSAR  vs  the  M  at  LASA,  for  the  events 

s  s 

listed  in  Appendix  I,  as  determined  from  the  LP  vertical  component  of  the  Rayleigh 
wave.  The  scatter  about  the  45°  line  in  Fig.  1  must  be  attributed  to  source -medium 
difference  between  the  sites,  observational  error,  local  effects  and  instrumental  dif¬ 
ference.  These  differences  must  account  for  an  M  at  LASA  0.  6  higher  than  Faldalen, 

s 

to  1. 1  lower  than  Faldalen.  Since  we  are  not  able  to  remove  the  marked  variation  be¬ 
tween  the  individual  events,  we  will  adopt  the  viewpoint  that  each  of  these  factors  con¬ 
tributes  to  one  of  two  types  of  differences.  One  averages  out  over  a  large  sample  of 
events;  the  other  introduces  a  bias  when  an  average  is  taken.  However,  the  populations 
of  events  which  we  use  are  not  large  enough  for  effects  with  zero  mean  to  completely 
average  out. 

Earthquakes  in  general  have  a  lobate  radiation  pattern.  The  observed  magnitude 
of  an  event  will  be  azimuthally  dependent  unless  a  correction  is  made.  Since  the  source 
type  and  fault  plane  dip  and  strike  are  not  known  for  these  events,  we  can  make  no  cor¬ 
rection  for  the  source  difference.  However,  when  the  differences  in  radiation  pattern 
seen  at  two  sites  are  averaged  over  a  large  number  of  events,  the  effect  tends  to  can¬ 
cel  out.  We  idealize  this  source  difference  as  having  zero  mean. 

The  travel  path  from  the  epicenter  of  each  event  to  Faldalen  or  LASA  may  be 
quite  different.  Surface  waves  can  be  scattered  by  continental  margins,  mountains,  and 
other  large  scale  changes  in  the  crust  and  upper  mantle  structure.  The  most  pronounced 
medium  effect  on  surface  waves  at  Norway  occurs  on  events  from  Asia.  Figure  2  dis¬ 
plays  a  striking  example  of  the  medium  effect  on  a  Rayleigh  wave  seen  at  Faldalen  from 
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the  Hindu  Kush  region.  Following  the  20-second  Rayleigh  wave  marked  by  the  cursors 
a  large  apparent  Airy  phase  arrives.  This  represents  energy  in  a  dispersive  wave 
train  which  has  not  been  dispersed.  Both  continental  and  oceanic  crust -upper  mantle 
models  have  Airy  phases.  However,  a  given  travel  path  crosses  varying  crustal  struc¬ 
tures.  Hence  the  period  at  which  an  Airy  phase  propagates  in  one  region  is  dispersed 
by  the  next.  The  presence  of  Airy  phases  from  Asian  events  indicates  that  the  travel 
path  from  these  events  encounters  very  nearly  identical  crust-upper  mantle  structure 
along  its  path  to  Faldalen.  Regions  with  events  exhibiting  Airy  phases  at  Faldalen  are 
E.  China  Sea,  Kurile  Islands  and  Hindu  Kush.  The  group  velocity  curve  of  Rayleigh 
waves  for  an  average  continental  structure  is  nearly  flat  at  20-second  period  indicating 
that  Airy  phases  occur  at  this  period.  The  group  velocity  curve  for  an  average  oceanic 
structure  is  quite  steep,  indicating  significant  dispersion  of  20-second  period  energy. 
Asian  events  have  purely  continental  paths  to  NORSAR  but  mixed  oceanic -continental 
paths  to  LASA.  Since  a  large  number  of  Asian  events  are  included  in  our  sample,  we 

expect  the  travel  path  effect  to  give  an  average  bias  of  larger  M  values  at  NORSAR 

s 

relative  to  LASA. 

The  local  effects  are  of  two  types :  differences  in  local  surface  wave  velocity, 
and  local  scattering  and  noise.  Though  a  surface  wave  at  two  sites  may  have  the  same 
energy,  the  amplitude  will  be  smaller  at  the  site  with  the  higher  Rayleigh  wave  velocity. 
The  local  scattering  of  events  and  interference  is  more  difficult  to  evaluate.  One  expects 
this  to  be  azimuthally  dependent.  The  local  effects  could  introduce  a  bias  into  the  aver¬ 
age  difference  in  Mg  between  two  sites. 

Hie  observational  error  we  assume  has  zero  mean.  Any  instrumental  differences 
would,  of  course,  bias  the  data.  Therefore,  any  mean  difference  between  LASA  and 
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Faldalen  must  be  explained  by  the  effect  of  differences  in  the  travel  path  to  each  site, 

local  effects  or  instrument  bias.  The  mean  difference  in  M  [M  (Faldalen)  —  M  (LASA)] 

s  s  s 

for  the  events  in  Fig.  1  is  0.  3.  The  solid  line  in  Fig.  1  is  shifted  by  this  mean  difference 
above  the  45°  line.  Far  more  travel  paths  to  Norway  than  to  LASA  cross  purely  conti¬ 
nental  structures.  Since  the  purely  continental  path  scatters  and  disperses  less  surface 
wave  energy,  we  expect  the  average  effect  of  travel  path  differences  will  give  larger  M 

s 

values  at  Faldalen.  For  this  reason  we  take  the  mean  difference  as  an  approximate  upper 
bound  on  the  instrumental  difference  and  the  local  effects. 

To  further  assess  the  extent  of  local  effects  and  compare  Mg  values  measured 
at  the  three  LP  sites  in  Norway  we  computed  the  mean  difference  in  Mg  for  those  13 
western  U.  S.  and  Mexican  events  (Appendix  II)  which  were  observed  at  more  than  one 
of  the  three  Norway  LP  sites  (Faldalen,  Trysil  and  Oyer).  Since  these  events  have 
nearly  the  same  travel  path,  only  instrument  bias  or  local  effects  will  plan  an  important 
role.  The  mean  difference  for  individual  sites  was: 

Oyer  —  Faldalen  Oyer  —  Trysil  Faldalen  —  Trysil 

0. 1  -  0.  3  -  0.  3 

[Note  that  (O  —  F)  +  (F  —  T)  /  (O  —  T),  since  all  sites  did  not  record  each  event.  ] 

We  find  differences  between  individual  sites  in  Norway  as  great  as  the  differences 
between  Faldalen  and  LASA.  These  differences  are  not  considered  unusual  and  are  prob¬ 
ably  due  to  azimuthally  dependent  local  effects  and  differences  in  local  crustal  structure, 
not  instrumental  differences.  An  approximate  upper  bound  for  a  mean  difference  ap¬ 
pears  to  be  0.  3  for  individual  sites.  One  must  be  aware,  however,  that  the  M  differ- 

s 

ence  between  two  sites  varies  erratically  from  event  to  event.  For  the  small  populations 
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of  events  used  in  this  study,  the  effects  we  assume  to  have  zero  mean  will  not  average 
out.  The  numerical  values  determined  for  intersite  mean  differences  are  only  approxi¬ 
mate.  Use  of  beamforming  techniques  using  the  entire  projected  LP  array  ought  to  par¬ 
tially  average  out  these  local  directional  differences  in  sensitivity  and  local  crustal 
structure  differences. 
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THE  M  VERSUS  m.  DISCRIMINANT 
 s_ b 

The  relationship  of  Mg  vsm^  as  determined  from  data  recorded  at  LASA  has 

4 

been  examined  for  four  regions  of  the  world  by  Capon,  et  al.  Since  this  discriminant 

has  been  successfully  applied  at  LASA,  we  want  to  compare  its  relative  effectiveness 

for  data  recorded  at  Norway.  Figure  3  shows  the  LASA  Mg  vs  m^  relationship  for  the 

four  regions  determined  by  Capon,  et  al.  Figure  4  presents  the  M  vs  determination 

at  Norway  for  26  events  of  various  depths  occurring  outside  the  western  U.  S.  and  Mexico. 

Unfortunately,  this  population  contains  only  one  presumed  explosion  and  only  an  upper 

bound  on  its  surface  wave  magnitude  was  obtained.  These  events  are  listed  in  Appendix  III. 

The  M^  vs  values  of  events  from  N.  E.  China,  Hindu  Kush  and  Turkey  lie 

near  the  explosion  region  in  Fig.  4.  Figure  5  presents  events  with  M  =  3.  0  and  depth 

s 

less  than  50  km  as  determined  by  the  U.  S.  Coast  and  Geodetic  Survey.  These  events 
are  listed  as  "shallow  events"  in  Appendix  III. 

The  exclusion  of  deep  events  and  undetermined  depth  events  from  the  M  vs  m. 

S  D 

plot  gives  a  significant  improvement  in  the  separation  of  the  earthquakes  and  the  one 
explosion.  Yet  for  this  preliminary  data  no  definite  conclusions  concerning  the  effective¬ 
ness  of  the  discriminant  Mg  vs  m^  at  Norway  can  be  drawn. 
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WESTERN  UNITED  STATES  EVENTS 


According  to  Liebermann  and  Pomeroy,  ^  the  Mg  —  values  determined  for 
western  U.  S.  events  are  anomalous  compared  with  the  world  average  for  both  explo¬ 
sions  and  earthquakes.  The  data  supporting  this  hypothesis  were  recorded  within  con¬ 
tinental  North  America  at  World  Wide  Standardized  Seismic  Stations  (WWSSS)  and 
Long  Range  Seismic  Measurements  (LRSM)  stations.  Their  separation  from  six 
earthquakes  and  five  Nevada  Test  Site  (NTS)  underground  explosions  was  complete. 

The  earthquakes  occurred  in  southern  California,  Nevada,  Montana  and  New  Mexico. 

To  confirm  that  similar  results  occur  for  measurements  at  teleseismic 
distances  outside  the  continental  boundaries  of  the  epicentral  region  we  have  analyzed 
eight  U.  S.  earthquakes,  three  Mexican  earthquakes,  four  NTS  explosions  and  one 
collapse  of  an  explosion  for  Mg  vs  m^.  The  m^  value  was  computed  using  the  Oyer 

subarray  and  the  M  value  from  a  long-period  beam  of  Oyer,  Trysil  and  Faldalen. 

s 

Figure  6  presents  the  Mg  vs  m^  plot.  The  separation  can  be  seen  to  be  complete  for 
this  set  of  data.  The  Mg  —  values  for  the  explosions  fall  within  the  world  wide 
range  of  earthquakes  and  the  earthquakes  of  the  western  U.  S.  appear  to  have  either 
an  Mg  value  above  the  world  average  or  a  low  m^  value.  The  dashed  line  would  separate 
earthquakes  and  explosions  for  the  data  of  Fig.  3  determined  at  LASA.  This  anomalous 
behavior  of  the  Mg  —  m^  value  for  western  U.  S.  events  clearly  confirms  indications 
from  Fig.  3  that  earthquakes  and  explosions  must  be  treated  on  a  regional  basis.  The 
Mg  discriminant  appears  to  discriminate  completely  for  shallow  events  (m^  >  5) 

which  have  been  regionalized. 
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CAUSE  OF  REGIONAL  VARIATIONS  IN  THE  M  vs  m,  DISCRIMINANT 

s  —  b 

The  necessity  that  the  Mg  vs  discriminant  be  applied  to  regionalized 
populations  of  events  could  severely  limit  its  use  as  a  discriminant.  Since  most  of 
the  earth  is  aseismic,  one  will  not  always  have  a  reference  population  of  earthquakes 
from  a  given  region  with  which  to  compare  future  events.  Clearly,  if  an  event  from 
such  a  region  has  an  Mg  vs  m^  value  similar  to  that  of  an  underground  explosion 
from  the  western  U.  S. ,  it  could  be  classified  as  an  explosion  using  the  western  U.  S. 
population  or  as  an  earthquake  using  the  Central  Asia  population.  The  problem  of  pre¬ 
dicting  a  priori  whether  a  given  region  is  similar  to  the  western  U.  S.  or  Central  Asia 
is  therefore  of  immediate  concern.  We  try  to  correlate  the  M  vs  m,  characteristics 
of  a  region  with  the  tectonic  character  and  the  upper  mantle  structure. 

The  cause  of  the  upward  shift  of  the  Mg  vs  values  for  events  from  the 
western  U.  S.  relative  to  Central  Asia  can  be  examined  from  two  viewpoints.  Either 
the  m^  value  from  this  region  is  anomalously  low  or  the  Mg  value  is  anomalously  high. 

Of  course  both  effects  could  be  acting.  We  will  now  examine  causes  that  can  realistically 
account  for  these  shifts. 

We  can  account  for  the  relative  attenuation  of  body  wave  originating  in  one 
region  relative  to  another  by  a  relative  difference  in  the  Q  of  the  low  velocity  layer. 

To  establish  an  upper  bound  on  the  effect  we  consider  a  low  velocity,  low  Q  region 
with  300  km  thickness.  In  Table  I  we  list  the  attenuation  experienced  by  waves  of 
three  frequencies  for  different  Q.  The  Q  effect  can  easily  account  for  the  shift  in  the 
Mg  vs  plot,  since  at  1  Hz  a  change  from  a  Q  of  200  to  50  in  this  layer  shifts 
downward  by  1.  5.  Such  differences  in  Q  may  very  well  exist  in  the  earth. 
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Q 


50 

100 

200 

500 

Frequency  0.  5 

.9  x  10"1 

.3 

.  6 

.8 

1.0 

.9  X  10"2 

.9  X10'1 

.3 

.6 

2.0 

.8  xio"3 

.9  x 10-2 

.9  x  10"1 

.4 

TABLE  I.  Attenuation  of  a  300  km -thick  layer 
(Transmitted  amplitude /incident  amplitude) 

As  a  preliminary  estimate  of  the  relative  effect  of  Q,  a  mean  spectrum  for  four 
populations  of  events  recorded  at  NORSAR  was  determined.  The  populations  consisted 
of  U.  S.  explosions,  USSR  presumed  explosions,  western  U.  S.  earthquakes  and 
worldwide  (excluding  the  U.  S. )  earthquakes.  The  spectrum  of  each  event  is  the  product 
of  the  source  spectrum,  the  response  of  the  source  region,  and  the  response  of  the 
receiver  region.  The  response  of  the  receiver  is  the  same  for  all  events,  while  the 
regional  difference  in  source  region  response  varies  between  the  populations.  By 
averaging  we  hope  to  obtain  an  average  source  spectra  for  the  earthquakes  and  for  the 
presumed  explosions,  though  this  may  not  be  possible.  Figure  7  shows  the  average 
spectrum  for  each  population.  The  higher  frequency  components  are  attenuated  more 
severely  on  the  U.  S.  spectra  than  the  others.  This  attenuation  would  occur  if  there  is 
a  regional  Q  difference.  However,  the  effect  we  are  observing  may  not  be  due  to  a 
difference  in  Q,  but  rather  regional  differences  in  source  spectra. 

The  relative  enhancement  of  the  surface  waves  in  one  region  relative  to  another 
is  attributed  to  a  different  effect.  Theory  predicts  that  an  explosive  source  in  an 
infinite  homogeneous  medium  will  not  generate  SH  waves.  Yet  it  is  well  known  that 
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nuclear  explosions  detonated  at  the  Nevada  Test  Site  (NTS)  often  generate  very  large 
Love  waves.  For  a  number  of  nuclear  explosions  Toksbz^  has  fitted  a  theoretical 
model  to  the  observed  data  using  amplitude  equalization  of  the  surface  waves.  By 
superposing  an  explosive  source  and  a  double  couple  source  one  can  fit  the  observed 
ratio  | L |  /  |r|  as  a  function  of  azimuth,  (where  |l|  is  the  maximum  amplitude  of 
the  Love  wave  and  |  R|  is  the  maximum  amplitude  of  the  Rayleigh  wave).  One  can  then 
determine  the  source  parameters  such  as  fault  plane  orientation  and  dip,  and  the  part 
of  the  source  due  to  double  couple  relative  to  the  explosive  component.  Toksdz  finds 
the  double  couple  component  contributing  as  much  as  three  times  the  explosive  com¬ 
ponent.  Similar  calculations  have  been  carried  out  by  the  author  for  later  NTS  shots, 
which  also  require  significant  double  couple  contributions.  The  source  mechanism 
of  an  underground  nuclear  explosion  near  NTS  must  allow  for  a  significant  earthquake 
component  as  well  as  an  explosive  component. 

The  source  of  the  SH  wave  generation  has  been  investigated  in  environments 

7 

where  source  conditions  can  be  determined  more  precisely.  Kisslinger  et  al  made 
measurements  near  the  source  of  chemical  explosions.  SH  waves  were  generated. 
When  surface  cracking  occurred,  it  had  a  preferred  orientation,  which  can  account 
for  the  asymmetrical  SH  radiation  pattern.  Later  model  experiments  reported  by 

g 

Kim  and  Kisslinger  give  more  conclusive  results.  Two-dimensional  prestressed 
models  were  excited  by  an  explosive  source;  shear  waves  were  generated  in  the  pre¬ 
stressed  models  but  absent  in  the  relaxed  ones.  Their  work  suggests'that  the  dominant 
role  for  S  wave  generation  is  attributable  to  the  strain  energy  release  rather  than  to 
directional  cracking.  "  The  S  wave  generation  is  increased  by  increasing  the  ambient 
stress  field  and  prefracturing  around  the  shot  hole. 
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If  one  assumes  the  models  behave  as  the  material  of  the  crust,  the  generation 

of  anomalously  high  surface  waves  by  underground  explosions  can  be  attributed  to 

9 

anomalously  high  ambient  stress  fields  in  the  upper  crust.  A  recent  article  synthe¬ 
sizes  data  from  many  areas  into  a  global  tectonic  picture.  The  model  distinguishes 
three  layered  regions  of  the  crust  and  mantle.  The  uppermost  mobile  region,  the 
lithosphere,  has  considerable  strength  on  a  geologic  time  scale.  The  lithosphere 
includes  generally  the  crust  and  upper  mantle  to  a  depth  of  100  km.  Below  the  litho¬ 
sphere  is  a  rather  weak  region,  the  asthenosphere ,  extending  to  several  hundred  kilo¬ 
meters  depths.  The  rest  of  the  mantle  is  the  mesosphere.  This  picture  includes  a 
continuation  of  an  oceanic  ridge  from  the  Gulf  of  California  along  the  San  Andreas 
Fault  to  the  west  coast  of  Oregon  and  Washington.  The  depth  of  seismic  activity  along 
the  San  Andreas  does  not  exceed  20  km  and  is  only  slightly  deeper  as  one  moves  away 
from  it.  Near  the  ridge  the  lithosphere  must  possess  a  large  ambient  stress  field  in 
the  upper  most  20  km  of  the  crust,  sufficient  to  push  apart  the  much  thicker  plates. 

An  explosion  detonated  near  the  ridge  of  this  lithospheric  block  will  release  substantial 
strain  energy. 

Earthquakes  occurring  in  the  western  U.  S.  must  also  release  more  strain 
energy  since  their  Mg  vs  plot  separates  from  that  of  the  explosions  in  the  region. 

To  examine  the  Mg  vs  m^  character  of  earthquakes  occurring  near  an  oceanic 
ridge,  the  Mg  vs_  m^  was  determined  for  a  population  of  39  earthquakes  near  the  mid- 
Atlantic  Ridge.  Figure  8  shows  this  plot.  All  but  two  of  the  events  lie  above  the 
Gutenbeig-Richter  line,  which  is  a  linear  fit  to  the  worldwide  average  Mg  v£  m^  for 
earthquakes.  These  results  support  the  belief  that  the  Mg  vs  m^  character  of  oceanic 
ridges  is  anomalously  high. 


12 


We  have  attributed  the  anomalous  M  vs  m,  behavior  of  the  western  U.  S.  to 

s  —  b 

two  possible  effects.  Further  investigation  is  necessary  to  determine  which  effect  is 
actually  contributing  the  most  to  the  shift  in  the  Mg  vs  m^  curve. 

To  determine  the  relative  Q  effect  of  the  western  U.  S.  with  other  regions  of 
the  world,  one  can  compare  the  spectra  of  a  deep  impulsive  event  recorded  at  many 
stations.  The  event  should  have  a  focal  depth  of  600  km  to  insure  that  the  first  arrival 
traverses  only  the  low  Q  region  at  the  various  receivers.  By  taking  records  from  both 
the  western  U.  S.  and  other  areas  one  could  evaluate  the  relative  attenuation  of  the 
high  frequencies.  Equivalently  one  has  a  measure  of  the  relative  Q  effect  between 
regions. 

The  relative  release  of  strain  energy  from  region  to  region  can  best  be  studied 
after  determining  the  Q  effect.  The  M  vs  m,  relationship  must  be  examined  for  a 
large  number  of  events  from  different  regions.  The  mean  focal  depth  for  these  regions 
should  be  determined  to  study  the  ambient  stress  distribution.  High  near-surface 
ambient  stress  field  should  correlate  well  with  high  Mg  vs  m^  values  for  events  from 
such  regions.  Large  stress  fields  at  depth  will  give  smaller  surface  waves  due  to  the 
lower  excitation  of  surface  waves  by  deep  sources.  If  the  release  of  strain  energy 
proves  to  be  the  cause  of  higher  Mg  vs  m^  values  for  nuclear  explosions  in  the  western 
U.  S.  ,  one  will  be  able  to  predict  a  priori  the  characteristic  M  vs  m,  values  for  a 
region  from  a  global  tectonic  model.  Furthermore,  one  will  have  the  assurance  that 
aseismic  regions  will  not  produce  anomalously  high  surface  waves  since  they  do  not 
possess  large  ambient  stress  fields. 
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EVENTS  USED  IN  COMPARISON  OF  M  AT  LASA  WITH  M  AT  FALDALEN 
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WESTERN  U.  S.  AND  MEXICAN  EVENTS 


Latitude  Longitude 

Shallow  Events 

Depth 

Region 
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16.  4  N 

92.1  E 

33 

Bay  of  Bengal 

4.3 

4.5 

49.8  N 

78.0  E 

0 
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<3.2 

5.5 

37.6  N 

114.9  E 

33 

N.  E.  China 

3.8 

5.7 

8.2  N 

38.2  W 

33 

Atlantic  Ridge 

4.  7 

5.7 

35.5  N 

22.5  E 

44 

South  of  Crete 

<3.3 

4.5 

46.4  N 

153.3  E 
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Kurile  Islands 
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132.2  E 
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FALDALEN 


Fig.  1.  Faldalen  Ms  V£  LASA  Ms  for  21  worldwide  earthquakes. 
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20-SEC  AIRY 

PERIOD  PHASE 

SURFACE 
WAVE 


Fig.  2.  An  Airy  phase  arriving  after  the  20 -second  period 
Rayleigh  wave  marked  by  the  cursor. 
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(a) 


Fig.  3.  Mg  vs^  plot  of  events  recorded  at  LASA  from  four  regions 
of  the  world  (a.  Aleutians,  b.  Central  Asia,  c.  Kuriles -Kamchatka, 
d.  Solomons -N.  Caledonia-N.  Hebrides). 
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SURFACE  WAVE  MAGNITUDE  (LASA) 


(b) 


Fig.  3.  Continued. 
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SURFACE  WAVE  MAGNITUDE  (LASA) 


(c) 


Fig.  3.  Continued. 
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SURFACE  WAVE  MAGNITUDE  (LASA) 


(d) 


Fig.  3.  Continued. 
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NORSAR 


~1 18-8-9  802~1 


SOLOMON 

ISLANDS 


□  PRESUMED  NUCLEAR  EXPLOSIONS 
•  SHALLOW  EARTHQUAKES 
O  EARTHQUAKES  WITH  CGS  DEPTH  >  50  km 


BAY  OF 
BENGAL 


TUNISIA 


CHILE-ARGENTINE 
BORDER 


N.  CHILE 
O 

E.  COAST 
KAMCHATKA 
O 


KURILES 

O  •  •  KURILES 
NORTH 
ATLANTI 
RIDGE 

SHIKOKU 


O  MINDANAO 


• 

CENTER 
r  MID-ATLANTIC 
RIDGE 


KERMADEC 

O 


. ANDAMAN 
ISLANDS 


N.E.  CHINA 


O 

HINDU  KUSH 


KIRGHIZ  O  O  BULGARIA 
S.  OF  CRETE 


YUGOSLAVIA' 

O 


E.  KAZAKH 


? 


o 

TURKEY 


YUGOSLAVIA 

O  O  ROMANIA 

I _ L 


4.5 

USCGS  m. 


5.0 


Fig.  4.  Mg  v^rrib  plot  of  25  earthquakes  and  one  presumed  explosion 
recorded  at  NORSAR. 
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NORSAR 


Fig.  5.  Ms  v£  mj,  plot  of  the  events  of  Fig.  4  whose  Ms  was  greater 
than  3.  0  and  depth  (as  determined  by  USCGS)  was  less  than  50  km. 
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Fig.  6.  Ms  V£  plot  of  nine  western  U.  S.  and  Mexican  earthquakes, 
one  NTS  collapse  and  four  NTS  explosions. 
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RELATIVE  AMPLITUDE 


Fig.  7.  Average  P  wave  spectra  recorded  at  NORSAR  from  U.  S.  and  non-U.  S. 
explosions  and  earthquakes.  The  spectrum  from  each  event  has  been  normalized 
by  the  square  root  of  the  sum  of  the  squares  of  the  spectral  components.  Each 
population  was  then  averaged  and  normalized  by  the  maximum  spectral  component 
for  comparison. 
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Fig.  8.  Mg  V£  plot  of  39  Mid-Atlantic  Ridge  events  recorded  at  LASA. 
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